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ABSTRACT 

We used archival Hubble Space Telescope WFC3 images to obtain the Luminosity 
Function of the remote globular cluster NGC 2419 from two magnitudes above the 
Horizontal Branch level down to ~ 3.0 magnitudes below the Turn Off point (to 
Mj ~ 6.4), approximately covering the range of initial stellar masses 0.5 Mq < m < 
0.9 Mq. The completeness-corrected Luminosity Function does not display any change 
of shape over the radial range covered by the WFC3 data, out to ~ 6 core radii (r c ), 
or, equivalently, to ~ 2 half-light radii. The Luminosity Function in this radial range 
is also identical to that obtained from ground based data at much larger distances 
from the cluster centre (12r c < R < 22r c ), in the magnitude range in which the two 
distributions overlap (Mj ^ 4.0). These results support the conclusion by Dalessandro 
et al. that there is no significant mass segregation among cluster stars, hence the stellar 
mass-to-light ratio remains constant with distance from the cluster centre. We fitted 
the observed Luminosity Function with theoretical counterparts with the proper age 
and metallicity from different sets of stellar evolution models and we consistently 
derive a total V band mass-to-light ratio 1.2 < M/Ly < 1.7, by extrapolating to the 
Hydrogen burning limit, with a best-fit value M/Ly = 1.5 ± 0.1. On the other hand, 
assuming that there are no cluster stars with m ^ 0.3 Mq, we establish a robust lower 
limit M/Ly > 0.8. These estimates provide useful constraints for dynamical models 
of the cluster that were forced to consider the stellar mass-to-light ratio as a (nearly) 
free parameter. 
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1 INTRODUCTION 

The remote and massive globular cluster (GC) NGC 2419 is 
a curious stellar system in several aspects. For example, it is 
the only bright (My — —9.5) Galactic GC lying at Galacto- 
centric distances larger than 30 kpc (see Fig. 8, in Galleti et 
aL]|2007 l. It has a half-light radius much larger than other 
GCs of the same luminosity ( Mackey & van den Bergh 2005 



et al.|2007 Brodie et al.|2011 I than to classical globulars. It 
displays a very extended bi-modal Horizontal Branch (HB 



Sandquist & Hess 2008 Dalessandro et al. 20081, possibly 



suggesting a significant spread in He abundance (di Cri 
|scienzo et al.|[2~011b[ ). A small spread in calcium and iron 
abundance has also been claimed, recently, based con Ca 



triplet spectroscopic surveys (Cohen et al. 2010 Ibata et 



Bellazzini|2007| and references therein) , more akin to stellar 



al.|2011a|). However the iron abundance of seven stars stud- 



nuclei or Ultra Compact Dwarf galaxies (UCD Evstigneva 



ied by Cohen, Huang & Kirby (20111 with high-resolution 



spectra appears remarkably homogeneous. Dalessan dro et| 
|al.| (|2008| D08 hereafter) used deep photometry over a field 



Based on observations made with the NASA/ESA Hubble 



Space Telescope, obtained at the Space Telescope Science Insti- 



tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555 
These observations are asso- ciated with program GO-11903 [P.I. 



covering the whole cluster to demonstrate that Blue Strag- 



gler stars (BSS, Bailyn 1995 1 in NGC 2419 have the same ra- 



J. Ka lirai] . 



T £j-mau: micnele.DellazzimaJoabo.mal.it 



dial distribution as Red Giant Branch (RGB) and HB stars 
(as well as the overall cluster light). The same behaviour 
has been observed only in the very peculiar massive cluster 



u) Cen ( Ferraro et al. 2006 I and in the sparse and remote 



M. Bellazzini et al. 



cluster Pal 14 (Beccari et al. 20111, while in all the other 



GCs, BSS are strongly centrally concentrated with respect 
to other stellar species, as a consequence of mass segregation 
driven by two-body relaxation (Ferraro & Lanzoni 2009). 



traced down to 0.10-0.15 Mq (i.e. near the Hydrogen burn- 
ing limit) with HST photometry reaching I ~ 25.0 (|Paust et 



In recent times, NGC 2419 has been the object of sev- 
eral kinematic studies aimed at deriving the mass distribu- 
tion within the cluster and at testing alternative theories of 
gravitation ( Baumgardt et al.|2009 Conroy, Loeb & Spergel 



2011 Ibata et al.||2011a|b l. This is due to several interest- 
ing (and in some case unique) properties of this system. Its 
large Galactocentric distance (Rgc — 95 kpc) implies that 
the effect of Galactic tidal forces on the cluster are negligi- 
ble, so that the motion of their stars at any distance from 
the cluster centre can be safely interpreted as driven only 



by the cluster potential and the law of gravitation ( Baum- 
|gardt et alT]|2009| [Ibata et al.||2011a| ). The total luminosity 
of NGC 2419 is ~ 25 times larger than any other Galactic 
GC lying in these remote regions. This means that the po- 
tential targets for estimating individual radial velocities, i.e. 
RGB stars reachable with medium resolution spectrographs 
on 8m class telescopes, are fairly abundant and can allow 
to sample adequately the velocity dispersion curve out to 
large distances (see Ibata et al.||2011a l, while in other dis- 
tant galactic GCs they are barely sufficient to derive a global 
value of the dispersion (see, e.g., Jordi et al.|2009 Sollima et 



aL]|2012 1. The gas- less nature and spherical symmetry typ- 



ical of old globulars greatly simplifies the modelling of the 
system. 

Since it is apparent that the cluster potential is domi- 
nated by the mass provided by stars, any independent con- 
straint on the stellar mass-to-light ratio (M/L) and its ra- 
dial behaviour can provide precious support to such studies. 
Since GCs stars are coeval and chemically homogen eous 
given a fixed mass-luminosity relation (MLR; see, e.g. Paust 



et al. 2009| and references therein), (a) the total stellar mass- 



to-light ratio depends only on the form of the present-day 
Mass Function (MItj except for the contribution of dark 
remnants, that is not dominant, see below), and (b) the ra- 
dial variation of the stellar M/L depends only on the degree 
of mass segregation within the cluster. Both issues can be ex- 
plored observationally by deriving a completeness-corrected 
Luminosity Function (LF) and looking for any radial varia- 
tion of the LF. This is the aim of the present analysis. We use 
the deepest images available for NGC 2419, obtained with 
the Wide Field Camera 3 (WFC3) on board of the Hubble 
Space Telescope (HST), to study the LF of the cluster in 
different radial ranges. 

The techniques to obtain reliable luminosity (and mass) 
functions in GCs, corrected for all the observational effects, 
are well established and have been widely used for several 
years (see |Paust et al.||2009| | 2010| |De Marchi et al.||2010 



|Bastian, Covey fe M eyer 2010[ for recent applications, re- 
view and references). While the MF of nearby clusters can be 



1 In the considered case any small i nhomogeneity j Cohen et al,| 
[20101 [Cohen, Huang fe Kirby|[2uTT| |di Criscienzo et al.||2011b| | 
should have negligible effects on the stellar mass-to-light ratio. 

2 Since we are interested in the current status of the cluster we 
consider only the present-day Mass Function, not the Initial Mass 
Function that is the subject of many studies in this field. In the 
following we refer to the present-day Mass Function simply as 
MF. 



al.|2010 l, the large distance to NGC 2419 limits the accessi- 
ble range of masses to in < 0.5 Mq even if our photometry 
reaches I ~ 27.5 (see below). Therefore the portion of the 
cluster MF providing the largest contribution to the stel- 
lar mass budgefj is out of reach of the present study. This 
will also limit our ability to detect radial variations, since 
a larger range of stellar masses sampled would increase the 
sensitivity to subtler effects. Moreover, the relatively small 
field covered by the WFC3 data also limits the explored ra- 
dial range to ~ 2 half-light radii (7%). Still, there are solid 
reasons to believe that the results that can be obtained from 
this analysis are worth the effort: 

(i) While the total V band mass-to-light ratio predicted 
by stellar evolution models are typically around ~ 2 (Krui- 
|jssen fc"M icskc 2009|, direct dynamical estimates as low as 
M/Ly = 0.4 — 0.5 are reported for some clusters (Mandu- 



shev et al.|1991||Pryor fc Meylan|1993||Strader et al.|2011 
In some cases this can be due to observational limitation! 
but, in general, it is attributed to preferential loss of low 
mass stars, due to the combined effects of mass-segregation 
and evaporation in a tidal field (sec Kruijsscn & Micskc 2009 
for discussion and references). Thus, even a best-effort esti- 
mate or a robust lower-limit based on the real MF would be 
valuable in narrowing the parameter space to be explored 
with dynamical models. 

(ii) Detailed theoretical analyses (e.g 



Baumgardt & 



Makino 2003 1 show that the effects of mass-segregation 



should be quite evident in the inner regions of clusters and 
the slope of the MF is observed to change also on mass 
ranges smaller than that considered here. 

(iii) We use the deep and wide-field photometry by D08 
to extend the analysis on the radial variation of the MF 
over the whole body of the cluster, albeit limited to an even 
smaller range of sampled mass. 

In the following we always adopt the cluster parameters 
listed in Table [I] that are the most reliable and up-to-date, 
from the literature. Since the MLR must be the same at 
all radii, comparing LFs in different radial ranges is strictly 
equivalent to comparing the MFs. Also, theoretical LFs can 
be produced for any assumed IMffjusing web tools that will 
be described below. Hence we avoid the transformation of 
the observed LFs into the corresponding MFs since it is not 
needed in the present context. Throughout the paper we 

3 To provide a quantitative idea, from an isochrone of age=12 
Gyr, [Fe/H ]=-2.1 and [a/Fe]+0.4 from the Dartmouth set | |Dotter| 

|et al.|2007) , a MS star at the Hydrogen burning limit has M/L ~ 
60, one with m=0.4 Mq has M/L ~ 11, and one at the Turn Off 
point M/L ~ 0.3, while a giant at the RGB bump has M/L ~ 
0.01 and at the RGB Tip M/L ~ 1.0004 (bolometric luminosity). 

4 For example, |Pryor fc Meylan| jl993| ) report M/Ly = 0.7 at 
the centre of NGC 2419 based on an estimate of the central ve- 
locity dispersion cro = 3.0 km s _1 , derived from a sample of 12 
stars presented by ( [Olszewski et al.|1993| ). However, recent stud- 
ies based on much larger samples, demonstrated that the actual 
value of cro is clearly higher than 6 km s~ x , implying M/Ly ~ 2 
I Baumgardt et al.|2009||Ibata et al.|20lia) . 

' in all these models the 1Mb' is then evolved into a LF/MF of 
the chosen age, thus allowing the comparison with the observed 
present-day LFs (see Renzini & Fusi Pecci|1988 I . 
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Table 1. Adopted cluster parameters. 



Parameter 



value 



Ref. 



(m — 
E(B 



M) 

-V) 



n-, 

Vt 

M v 

[Fe/H] 

[a/Fe] 

RF606Wr\ 

Rfshw 



19.71±0.08 


0.08±0.01 


27.6" 


56.3" 


10.47±0.07 


-2.1±0.1 


+0.^ 


2.874 


1.894 



di Criscicnzo ct al. (2011a) 



di (Jris cicnzo ct al. (2011a) 
2Ulla| 



Ibata ct al 



Ibata et al. 



Bella 



2011a 



(2007) 



Cohen, Huang fc Kirby 1 2011 1 



Cohen, Huang fc K irby (2011) 
Marigo et al.| (2OO8 171 



Marigo et al 



(2008) 



a For their best-fit model (#17) 
b 



by 
c 1 



In excellent agreement with the indepedent estimate 
BellazzTni| | |2u07[ |. 
'torn the values listed above. 

d Average of [O/Fe], [Mg/Fe], [Si/Fe], and [Ca/Fe]. 

e Wher e R x = A X /E(B - V). 

J From http://stev.oapd.inaf.it/cgi-bin/cmd for 



the WFC3 VEGAMAG system 

will consider LFs in the F814W VEGAMAG HST-WFC3 
passband. In addition to the lower sensitivity to interstellar 
reddening with respect to F606W, we found that F814W 
can be safely transformed into Cousins' I, thus allowing 
an easy comparison with ground-based data. In Fig. IT] we 
compare F606W, F814W magnitudes with Johnson-Kron- 
Cousins V,I, as a function of colour, for more than 150 stars 
in common between our WFC3 photometry of NGC 2419 



and the set of secondary standards by Stetson (2000 120051). 
While the transformation between V and F606W display a 
strong (and possibly non linear) dependence on colour and 
a large scatter at any colour, I and F814W appear identical 
independently of colour, within the uncertainties. For this 
reason we will use F814W and I interchangeably in the fol- 
lowing. Note that the comparison of Fig.lllis performed with 
stars lying within 120" of the cluster centre, where crowding 
is expected to provide significant contribution to the r.m.s. 
scatter. Therefore, it is likely that magnitudes in the two 
filters do coincide to within less than ~ 0.02 mag, at least 
for metal-poor stars in this range of colour. 

The paper is organised as follows. In Sect. 2 we present 
the observational material and describe the data reduction, 
including the artificial stars experiments. In Sect. 3 we will 
derive the Luminosity Function at different radii and we will 
show that it is the same everywhere. Sect. 4 is devoted to 
obtain an estimate of the stellar mass-to-light ratio by means 
of comparison with different theoretical models. The main 
results of the analysis are summarised in Sect. 5. 



2 OBSERVATIONS AND DATA REDUCTIONS 

The data-set used for the present analysis consists of images 
obtained with the UVIS channel of the Wide Field Camera 3 
(WFC3) aboard the HST. These images are part of a photo- 
metric calibration program (Prop. ID: 11903; PI: J. Kalirai) 
so in principle they are not intended for LFs studies. Never- 
theless thanks to the way in which they have been planned 
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Figure 1. Transformations between the WFC3 VEGAMAG 
F606W and F814W magnitudes and the Johnson-Kron-Cousins 
V and I magnitudes, obtained from stars in common between our 
WFC3 sample and the secondary photometric standard sample 
by |Stetson| ( j2000| ) . Stars marked with a cross have been rejected, 
since they are outside the ±0.1 mag band (dotted lines) around 
the best fit (44 rejected stars in V-F606W and 31 in I-F814W). 



and the extraordinary efficiency of the camera, they result 
to be the deepest images available for NGC 2419. 
Given the aim of the present analysis we used only F606W 
and F814W images which allow one to reach fainter mag- 
nitudes than images with similar exposure times taken in 
different passbands and a straightforward comparison with 
theoretical models. The data-set used is thus composed by 
two images for each filter with the same exposure time 
t exp = 400 s for F 606W and t exp = 650 s for F814W. 
Images are perfectly aligned, with the cluster centre lying in 
chip #1, at about 10" from the gap which splits the camera 
in two twin chips. 

The photometric analysis has been performed on calibrated 
(i.e, bias and dark subtracted, and flat-field corrected) sin- 
gle exposure images (-fit) processed by the standard HST 
pipeline calwf3. To these images we applied the most up- 
dated Pixel Area Map (PAM) corrections by using standard 
IRAF tasks. 

In order to obtain a reference frame cleaned of cosmic-rays 
and with an enhanced signal to noise ratio, we combined all 
the images in our data-set with the IRAF task imcombine 
applying a cosmic-ray rejection algorithm. 
For each filter and chip, a large number of relatively bright 
and isolated stars have been selected in order to model prop- 
erly the Point Spread Function (PSF) all over the Field of 



View. The PSF was modeled with a Moffat ( 1969 1 function 



the parameters of which (a, /3) have been allowed to vary 
linearly as a function of position. 



We used DAOPHOTII/ALLSTAR ( |Stetson||1987| ) to 
search for sources at 2<r above the background on the 
cleaned reference frame, thus obtaining a master list of 
bona fide stars. Then, at the corresponding positions of 
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Figure 2. Colour Magnitude Diagrams of NGC 2419 from HST- 
WFC3 data, in the VEGAMAG system. Left panel: observed 
CMD. Right panel: CMD corrected for distance and reddening. 
The box adopted to select RGB+MS stars and HB stars arc plot- 
ted as continuos lines. 



these stars in the single images, a fit was attempted with 
DAOPHOTII/ALLFRAME. For the stars recovered at 
the end of the procedure and for each filter the resulting 
magnitudes were averaged and centroids positions combined 
by using DAOMATCH and DAOMASTER. 

The instrumental magnitudes thus obtained were trans- 
formed into the VEGAMAG photometric system using pre- 
scriptions and photometric zeropoints reported in the HST 
web-sitej_] 

In order to assign absolute coordinates to each star of our 
final catalog, we corrected for the heavy geometric distor- 
tions which affects the WFC3, by applying the equations by 
Bellini & Bedin|(|2009|. We used the stars in the data-set by 
D08 as secondary astrometric standards. We used CataX- 
corqjto transform the WFC3 instrumental coordinates into 
Equatorial coordinate systems by using a few thousand stars 
in common with D08. 

In the left panel of Fig. [2] we present the Color Mag- 
nitude Diagram (CMD) for the whole sample. The cluster 
Main Sequence dominates the diagram, with a Turn Off 
(TO) point around F814W ~ 23.5 and reaching the limiting 
magnitude of the photometry at F814W ~ 27.5. A sparse 
plume of BSS is clearly visible around F606W - F814W ~ 
0.0 and F814W ~ 22.5. Stars on the steep Red Giant 
Branch, typical of metal-poor clusters, are reliably measured 
up to F814W ~ 17.5, where saturation sets in. The ex- 
tended and bimodal HB reaches F8UW ~ 25.5 (see D08, 



6 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 

7 CataXcorr is a code aimed at cross-correlating catalogues and 
finding astrometric solutions, developed by P. Montegriffo at 
INAF - Osservatorio Astronomico di Bologna, and successfully 
used by our group for the past 10 years. 



di Criscienzo et al.|2011b for details and discussion). In the 
right panel of Fig. |2j magnitudes have been corrected for 
distance and reddening, and we show the selection boxes 
that we adopted to select likely cluster members for the 
RGB+MS and HB+AGB (Asymptotic Giant Branch). We 
use RGB+MS stars to derive the cluster LF in Sect. 3, while 
the HB+AGB sample is used in Sect. 4.1 to estimate the im- 
pact of these stars on the global mass-to-light ratio. 

All over the paper we will consider the sample of stars 
having the DAOPHOT sharpness parameter in the range 
—0.3 < SHARP < 0.3; this selection criteria is found to 
clean the sample from sources whose shape is clearly dif- 
ferent from the bulk of genuine stars without significantly 
affecting the overall sample size. Among the 156281 stars 
lying within the RGB+MS selection box, 139209 fulfil the 
-0.3 < SHARP < 0.3 criterion. 

The contamination by foreground Galactic stars is ex- 
pected to have little impact over such a small field of view 
as considered here (162" x 162" ~ 0.002 deg 2 ). The pre- 



dictions of the TRILEGAL Galactic model ( Girardi et al 



2005 1 and our analysis of the stellar population in the 
external region (beyond R ~ 600") of the D08 field in- 
dicates that less than ten Galactic stars are expected to 
contaminate our sample, all of them being brighter than 
F814W ~ 23.4 (M FS i4w ^ 3.5). In particular, from the 
D08 dataset, we estimate that the fraction of Galactic stars 
in the sample selected for the analysis of the LF in the range 
17.9 < F814W < 23.4 (-2.0 < M F8 nw < 3.5) is less than 
0.01 per cent. At fainter magnitudes unresolved background 
galaxies are, by far, the main source of contamination. Form 
the D08 dataset we estimate that the contamination by 
unresolved galaxies in the range 23.4 < F814W < 25.4 
(3.5 < Mf814w ^ 5.5) is about 0.2 per cent. This may 
slightly under-estimate the actual contamination because 
it is not corrected for completeness, since we do not have 
such correction for the D08 sample; on the other hand it 
may overestimate the contamination because many of these 
sources that appear star-like in the Subaru images used by 
D08 would be recognised as non-stellar in the WFC3 images. 
The limiting magnitude of the D08 sample (see Fig. FFl be- 
low) prevented us to study the contamination fainter than 
F814W ~ 25.4. From th e publicly available ca talogue of the 
HST Ultra Deep FielcQ ( |Beckwith et al.|2005[ ) we estimated 
that, for F814W > 25.4 to the limiting magnitude of our 
photometry, the contamination by background galaxies is 
lower than a few per cent at any magnitude, i.e. smaller than 
the observational uncertainties in all the LFs we present in 
the following analysis. 



2.1 Artificial Stars Experiments 

A crucial step in determining an accurate LF is a proper 
estimation of photometric completeness. This aspect is par- 
ticularly important in the case of NGC 2419, which appears 
to be a remarkably crowded cluster because of its distance. 
Artificial stars experiments were performed following 



the method described by Bellazzini et al. (2002) and Bec- 
cari et al. (20101. We generated a catalog of simulated stars 



with a F8l4Wi n magnitude extracted from a LF modelled 



http : //archive . stsci . edu/pub/hlsp/udf /acs-wf c/h_udf _wf c_Vl_cat . txt 
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Figure 3. Colour Magnitude Diagrams and Luminosity Func- 
tions for artificial stars. Left panel: CMD from input magnitudes. 
Middle Panel: CMD for output magnitudes. Right panel: Lumi- 
nosity Function for the input (continuous histogram) and output 
(dashed histogram) samples. 



to reproduce the observed LF in that band and extrapo- 



lated beyond the limiting magnitude (see |Bellazzini et al.| 
2002 ). Then to each star extracted from the LF we assigned 
a F&QQWin magnitude by means of an interpolation along 
the mean ridge line of the cluster (left panel of Figure 13k. 

Artificial stars were added to real images (including 
the reference frame) by using the DAOPHOTII/ADDSTAR 
software. In order to avoid "artificial crowding", stars were 
placed into the images following a regular grid denned by 
15 pixels cells (6-7 times the usual FWHM of stars for these 
images) in which only one artificial star for each run was 
allowed to lie. 

The photometric reduction process used for the artifi- 
cial stars experiment is exactly the same as described in Sec- 
tion^ Those stars recovered after the photometric analysis 
have also a F606W O ut and F8UW ut- More than 200000 
stars have been simulated for each chip. The photometric 
completeness C/ is defined as the ratio between the number 
of stars recovered at the end of the procedure and the total 
number of stars actually simulated. 

The input and output CMDs and LFs of the artificial 
stars sample are compared in Fig. [3] Note the close similar- 
ity of the synthetic CMD with real one plotted in Fig. [2] 
Throughout the paper we always adopt the same selection 
criteria for observed and artificial stars. 

In Figureplwe compare the residuals between input and 
output magnitudes for both filters. As expected, the distri- 
bution is not symmetrical: a significant number of stars have 
been recovered with an output magnitude brighter than the 
assigned input of a quantity larger of their typical photomet- 
ric error. This effect is due to those stars that blended with 
real stars with similar (or larger) luminosity. This effect has 
been taken into account properly when building the LF (see 
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Figure 4. Distribution of the input-output magnitude differ- 
ences, as a function of output Mpsi4W magnitude. 



2.2 Completeness Functions 

In the following we will consider the sample of stars en- 
closed in the MS+RGB selection box displayed in Fig. [2] 
and the analysis will be performed using absolute magni- 
tudes (Mf814w)- 

In Fig. [5] we plot the completeness factors as a function 
of magnitude for several ranges of distance from the clus- 
ter centre (R). The crucial conclusion that can be drawn 
from this plot is that at any given magnitude fainter than 
Mf8uw = 2.0, the degree of completeness is subject to sig- 
nificant variations on scales as small as ~ 10". Conversely, 
for Mfsuw ^ 2.0, the completeness is 100% at any radius; 
this is true also for Mfsuw ^ 3.0, at least for R > 20.0". 
Usually, in nearby GCs, the completeness is nearly constant 
over relatively large radial intervals, so that a unique correc- 
tion can be applied to the LF from stars in that range (see, 
e.g. Bellazzini et al. 2002). This is clearly not the case for 



the remote cluster considered in this study. To account for 
the fast variation of Cj (m) with radius and for the effects of 
bin migration along the LF due to blending convolved with 
photometric error (see Fig.pl, we adopted the following cor- 
rection procedure acting on a star-by-star basis, inspired by 



Bailyn et al. (19911 and |Bellazzini et al. (20021 



Section 2.2 \ 



(i) For each star with Mfsi4w > 2.0, and having 
Mfsuw = m o an d R = ro, we derive the Complete- 
ness as a function of magnitude, with the associated un- 
certainties, from the sample of artificial stars lying within 
ro — 5.0" < R < ro + 5.0", with the same binning as in 
Fig-0 

(ii) We interpolate on the derived local C/(mrj with 



9 That would be more properly denoted as Cf(m,r), to empha- 
sise its local nature. However we prefer to drop the r argument 
to avoid a too heavy notation. 
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splines to obtain the completeness and the associated uncer- 
tainty at mo, C/(mo) and errCf(mo). The construction of 
the whole Cf (m) allows a more reliable (less noisy) estimate 
of Cf(mo) than extracting only the stars in the magnitude 
bin of the considered star, since the overall shape of the 
Cf(m) contributes to the derivation of the actual C/(mo) 
value. 

(iii) We extract at random a Gaussian deviate Gdev 
and we assign to the star the completeness factor C/> = 
Cf(mo) + Gdev x errCf(mo), with the additional constraints 
that if the value of C/,+ exceeds 1.0 its is set equal to 1.0, and 
if it is lower than 0.0 it is set to 0.0, so that 0.0 Sj C/,* ^ 1.0. 
This step takes into account the effect of the uncertainty in 
the determination of C/(mo). 

(iv) The subsample of stars having mo — 0.5 < m out < 
mo +0.5 is extracted from the sample of artificial stars lying 
within r — 5.0" < R < r + 5.0". We compute Am = 
m out — rrii n for these stars and we chose one of them at 
random. Its Am is then assigned to the considered real star 
(Am*). The magnitude of the real star is then corrected by 
the effect of bin-migration (blending + photometric error) 
by subtracting Am* to its observed magnitude, i.e. m* = 
mo — Am*. While the correction it is clearly unjustified on an 
individual star basis, on a large sample as the one considered 
here, it should provide a sound statistical correction. 

Once the above four steps are performed for all the stars 
in the sample, each observed star has a migration-corrected 
magnitude m* and an individual completeness value C/,*. 
Both of them have been derived from stars lying in the same 
radial and magnitude range of the considered star, so that 
the variations of these effects as a function of m and R are 
properly accounted for. The only limitation is due to the 
adopted size of the radius and magnitude bins, that should 
be sufficiently large to ensure robust determination of m* 
and C/,*. 

(v) Now the production of the corrected LF for any radial 
range is straightforward. One has simply to divide the stars 
lying in that range in magnitude bins, according to their 
m*, and then, instead of counting the stars falling within 
each bin i, the inverse of t heir completeness fa ctors must be 
summed Nj = y^ ^— (see 
we also record the average C~j 
only bins having (C/) > 0.5. 

(vi) To minimise the shot noise introduced by the random 
extractions included in the procedure we produced ten reali- 
sations of the catalog with associated m* and C/,*. Each LF 
presented below is obtained by averaging each Ni over the 
ten realisations. The associated standard deviation is added 
in quadrature to the Poisson noise of each bin to obtain a 
final uncertainty on (Ni). 



3 THE LUMINOSITY FUNCTION 

In Fig. Rj] we show the corrected LFs in four radial ranges 
sampled by our WFC3 data. These approximately corre- 
spond to Q.0-0.5rh, 0.5-1. Or^, 1.0-1.5rh, and 1.5-2. Or/,, from 
top to bottom, respectively. The faint branch of the LFs 
(Mfsi4w > 2.0), plotted as open circles, have been derived 
as described above, while for the bright branch (Mfsuw < 
3.0, open squares) we produced the LFs as ordinary his- 



Bailyn et al.ll991 1. For each bin 



and we retain for the final LF 




M„,„ 



Figure 5. Completeness as a function of magnitude for several 
radial ranges (the leftmost curve is for 0" < R. ^ 10", the right- 
most is for 80" < R ^ 120". The lines connecting the points 
are spline interpolations through the data. The lines through the 
edges of the error-bars (reported for the completeness functions 
of four radial ranges, as examples) show how the uncertainties 
can be reliably interpolated with splines. 



tograms since, as mentioned above, in the Mpsi4w < 2.0 
range Cf = 1.0 at all radii. We extended the bright branch 
down to Mfsi4,w < 3.0 (Cf = 1.0 anywhere except for the 
innermost regions, where it is > 90%) to verify that the faint 
and bright branches are in agreement in the region of over- 
lap. This is indeed the case without any adjustment of the 
normalisation, thus demonstrating that the adopted correc- 
tion procedure gives self-consistent results. 

The LFs reaches different magnitude linrits on the faint 
end. Due to increasing crowding toward the centre, the 
(Cf) — 0.5 threshold occurs at brighter magnitudes in 
the inner regions. The LF in the range 60" < R < 90" 
reaches the same limit as the outermost LF (Mfsuw = 6.2 
at the bin centre, implying the inclusion of stars down to 
Mfsuw = 6.4, that is the faintest limit of our LFs) but has 
been obtained from a much larger sample. For this reason 
we take it as a reference for comparison with LFs in other 
radial ranges and we list it in Tab. [3] The line connecting 
the points of the reference LF is plotted as a dotted line 
in the panels showing the LF in the other radial ranges in 
Fig. [51 after normalisation to the four points in the range 
0.5 < Mfsuw < 2.0. This approach to nornralisation is 
applied in all cases of comparison between empirical or the- 
oretical LFs in the following analysis. 

Our LF samples the cluster population from ~ 2 mag 
above the HB to 2-3 mag below the Main Sequence TO, 
depending on the considered radial range. Adopting an 
isochrone from the Dartmouth set ( jDotter et al.|2 007), with 
[Fe/H]=-2.1, [a/Fe]=+0.4, Y=0.2455 and age = 12 Gyr 
(that provides a good fit to the observed CMD, see Fig.pl be- 
low), the magnitude range covered by the reference LF cor- 
responds to the mass range between 0.81 Mq and 0.53 Mq. 
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Table 2. Luminosity Function of NGC 2419 in the 
range 60" < R < 90". 



A^814W 


logisg 


err (log N) 


brand 


b 


-1.3 


1.431 


0.083 


1 


-0.9 


1.477 


0.079 


1 


-0.5 


1.362 


0.090 


1 


-0.1 


1.519 


0.076 


1 


0.3 


1.716 


0.060 


1 


0.7 


1.869 


0.050 


1 


1.1 


1.973 


0.045 


1 


1.5 


2.017 


0.043 


1 


1.9 


2.260 


0.032 


1 


2.2 


2.331 


0.039 


2 


2.6 


2.574 


0.032 


2 


3.0 


2.882 


0.028 


2 


3.4 


3.163 


0.026 


2 


3.8 


3.314 


0.025 


2 


4.2 


3.458 


0.025 


2 


4.6 


3.571 


0.025 


2 


5.0 


3.670 


0.025 


2 


5.4 


3.746 


0.026 


2 


5.8 


3.830 


0.029 


2 


6.2 


3.966 


0.038 


2 





a Computed on 0.4 mag bins. 

" 1 = bright branch; 2 = faint branch. 



The key result of Fig. [6] is that, in the range of magni- 
tudes/masses covered by our LFs, there is no difference in 
the LFs obtained in different radial ranges, over O.Orj, < R < 
2. Or/,, or, equivalently, 0.0r c < R < 6.0r c . The shape of var- 
ious LFs is clearly the same and any subtle difference lies 
well within the combined uncertainty of the LFs and of the 
adopted normalisation. 

The comparison presented in Fig. [5] suffers from two 
main limitations. The first is that a large range of stellar 
masses is beyond the reach of our LFs, i.e. from 0.5 Mq to 
the hydrogen burning limit at ~ 0.1 Mq. Given the large 
distance of NGC 2419 this would be very difficult to over- 
come. A simple experiment with the WFC3 Exposure Time 
Calculator (ETC) shows that to detect stars two magnitudes 
fainter than our limit with S/N=3, more than 70 h of ex- 
posure time per filter are required. Moreover, pushing the 
limit of the LF two mag fainter than ours would still imply a 
mass limit of 0.3 Mq, still significantly above the H burning 
limit. Hence, the extension of the LF/MF to a significantly 
larger mass range seems impractical, at the present time. 
The second limitation arises from the radial range covered 
by the WFC3 data, since the cluster s extends at least out to 
7.6 r h | |Bellazzini|2007| P~ 



426" 



This can be partially 



n 

overcome by using the deep wide-field ground-based pho- 
tometry obtained by D08 from SuprimeCam@Subaru data, 
taking advantage of the equivalence between F814W and 
I photometry that we discussed in Sect. 1. In the upper 
panels of Fig. W\ we display the adopted selections for this 



10 The best-fit model by |lbata et al.| J2011a] has r t = 781", 
however, beyond R = 410" the Surface Brightness profile fall 
below fiy — 29.0 mag/arcsec 2 , implying a surface density of stars 
that is probably too low to derive statistically sound LFs with 
HST instruments. 
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Figure 6. Luminosity Function of NGC 2419 for four radial 
ranges (note that R = 60" corresponds approximately to the 
half-light radius of the cluster). The faint part of the LF (open 
circles) have been obtained applying the full Completeness Cor- 
rection, the bright part (open squares) refers to the magnitude 
range where Cf = 1.0 at any radius. The two faintest points of 
the bright LF (x symbols) overlap with the faint LF, showing 
that the two branches of the LF are fully compatible, i.e. there is 
no mismatch in the normalisation. The continuous line connecting 
the data points of the 60" < R ^ 90" is compared (dotted lines) 
to the LFs in all the other panels, after normalisation to the four 
points in the range 0.5 < Mpsuw < 2.0. The thin dotted lines 
enclose the ±2<r range of the uncertainty on the normalisation 
factor. 



sample. Since in the regime of moderate crowding the com- 
pleteness typically reaches its flat branch about 2 magni- 
tudes above the limiting magnitude of the sample, we con- 
sidered only stars lying above this limit (Mi < 4.0) and 
we derived the LF without any completeness correction. We 
focussed on the radial range 240" < R ^ 426" correspond- 
ing to 4.3rh < R < 7.6rh, i.e. sufficiently external to avoid 
the most crowded region and sufficiently large to obtain a 
reliable LF of the cluster outskirts. 

In the lower panel of Fig. [7] we compare the LF from the 
Subaru data to our reference LF (dotted line), with the usual 
normalisation. Also in this case the LFs in these two widely 
different radial ranges clearly display the same shape. We 
verified that contamination from foreground sources and/or 
background unresolved galaxies does not significantly af- 
fect this comparison (see Sect. 2). Hence, our analysis of 
the radial behaviour of the LF fully supports the conclu- 
sion reached by D08 using the BSS as tracer: the stars of 
NGC 2419 show no sign of mass-segregation all over the 
whole body of the cluster (see Baumgardt et aL]|2009 Da- 
|lessandro et al.|2008| for discussion). This evidence strongly 
supports the hypothesis of constant M/L with radius that 
has been adopted by Baumgardt et al. ( 2009 1 and Ibata et 
|al.| ( |2011a|b| ) in their analysis of the cluster dynamics. 

However, there is a caveat that must be mentioned. 



Our conclusion is limited to light-emitting stars: dark rem- 
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Figure 7. Upper left panel: map of the combined ACS + 
SUBARU-Suprimecam sample by D08 (black dots; only stars with 
Mj ^ 4.0 have been plotted) with the WFC3 sample superposed 
(light gray area). The circles are centred on the centre of the 
cluster and have radius=60", 120", 240", and 426". Upper right 
panel: CMD from the same sample with the box adopted to se- 
lect cluster stars for the LF. The long dashed lines mark the 
limiting magnitude and Mi = 4.0, two magnitudes above that 
level. Lower panel: LF for cluster stars with Mf < 4.0 in the ra- 
dial range corresponding to 4.3r^ < R < 7.6r^. The dotted line 
is the LF from the WFC3 sample in the range 60" < R < 90" 
(see Fig. [61, after normalisation to the four points in the range 
0.5 < Mpsi4W < 2.0. The thin dotted lines enclose the ±2<j 
range of the uncertainty on the normalisation factor. 



nants much more massive than the average cluster star (like 
black holes - BH - and neutron stars - NS) may be cen- 
trally segregated ( |Spitzer|1969| ), i.e. they may have a radial 
distribution different from the stars. The possibility that 
this would have a significant effect on the overall mass dis- 
tribution does not seems likely but it cannot be excluded. 
Given the present day status of the cluster and using the 



Equation 34 by Mouri & Taniguchi (20021, we find that the 



mass-segregation timescale (r ms ) in NGC 2419 is as large as 
11 Gyr for m — 1.4 Mq NS, and it becomes significantly 
smaller than the cluster age only for BHs with m > 10 Mq 
(j ms — 1.6 Gyr). Such high mass remnants should be quite 
rare and are not expected to provide large contribution to 
the overall mass budget. In Sect. [4~l] and Sect. [5] below we 
will discuss in more detail the contribution of dark remnants 
(including white dwarfs) to the total cluster mass. 



4 COMPARISON WITH THEORETICAL 
MODELS 

To constrain the global mass-to-light ratio we adopt the fol- 
lowing approach: 

• We fit the observed LFs with theoretical counterparts 
from evolutionary models with appropriate physical param- 
eters. There are now dedicated web tools that allow one to 



Figure 8. Comparison of the observed CMD with theoretical 
isochrones from the set by |Dotter et al.| | |2007[ l with [Fe/H] =- 
2.1. Left panel: isochrones with [a/Fe]=+0.2, and age=12.0, 13.0, 
and 14.0 Gyr, from left to right. Right panel: isochrones with 
[a/Fe]=+0.4, and age=11.0, 12.0, and 13.0 Gyr, from left to right. 



produce LFs of the desired form (e.g. choosing the exponent 
x of the power law IMF dN oc m x dm; where x = —2.35 
corresponds to the popular IMF of Salpctcr 1955) from the- 
oretical isochrones of given chemical composition and age, 
so it is easy to assemble a grid of theoretical LFs to compare 
with observations. 

• The best fit theoretical LF can be numerically inte- 
grated to obtain the total mass and the total luminosity (in 
any given passband), thus obtaining the stellar M/L inde- 
pendent of the normalisation. The effect of dark remnants 
can be taken into account with suitable assumptions (see 
|Kruijssen|2"009l ). 



al 



We adopt the Dartmouth theoretical models (Dotter et 
2007|P^ as our reference set. There are several reasons for 
this choice: (a) theoretical LFs and isochrones are provided 
down to the H-burning limit, in the WFC3 system, (b) the 
desired LFs can be produced with a web toop~land (c) this 



is the same set used by Paust et al. (2009 20101 in their 



exhaustive study of GCs mass functions, hence we have an 
homogeneous study to compare with. In the present context 



it is relevant to recall that Paust et al. ( 2009 1 showed that 



different sets of theoretical models provide MLRs for metal- 
poor stars in good agreement with each other, hence this is 
not expected to be a significant source of uncertainty in the 
determination of M/L ratios. 

In Fig. [8] we compare theoretical isochrones from the 
Dartmouth set with the observed CMD of NGC 2419. All 
isochrones have [Fe/H] = —2.1, to match the spectroscopic 



estimate by Cohen, Huang & Kirby (20111. Those in the 



left panel have [a/Fe] = +0.2, while those on the right have 



http : //stellar . dartmouth . edu/~models/ 



http : //stellar . dartmouth . edu/~models/isolf . html 
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[a/Fe] — +0.4. The latter value is equal to the measured av- 
erage reported in Table [I] however, since the average abun- 
dance of each individual a element can be quite different 
from this (even if always > 0.0 Cohen, Huang & Kirby 2011) 
it seems worth to explore the effect of different assumptions 
on [a/Fe] . It is quite clear that satisfactory fits of the MS 
and TO regions of the CMD can be obtained with both 
sets of isochrones. The well known age-metallicity degener- 
acy manifests in the fact that a slightly younger isochrone 
is required to fit the CMD with the set with the larger 
overall metallicity (i.e. the one with [a/Fe] = +0.4), 12 
Gyr instead of 13 Gyr. We have verified that the choice be- 
tween the two sets does not have any significant impact on 
the LF fit and the consequent M/L estimate. In the range 
of age and global metallicity [M/H] rH that is relevant for 
the present case the effects of these parameters on the LF 
shape are quite weak, the dominant parameter is the expo- 
nent of the MF power law. In the following we will adopt 
[Fe/H] = -2.1, [a/Fe] = +0.4 and age=12.0 Gyr as the 
standard set of values for the Dartmouth LFs to be com- 
pared with observations. When using other sets of theoreti- 
cal models, the BASTI set flPietrinferni et al.|2004| P~ 



(14 



or the 



Padua set (Marigo et al.|2008fl la we keep the global metal- 



licity fixed and we slightly adjust the age to fit the observed 
LFs. We recall that the normalisation between observed and 
theoretical LFs is always performed on the four points in the 
range 0.5 < M FS1AW < 2.0. 



Paust et al. (20101 conclude that the global present- 



day MF for all the GCs in their sample is well matched 
by a single power-law, in the range of masses between ~ 
0.8 Mq to 0.2 - 0.3 Mq. While it is like ly that a turnover 



in the MF occurs around these masses (De Marchi et al 



2010 Bastian, Covey & Meyer 2010 1, our data barely reaches 



m = 0.5 Mq, hence a comparison with single power-law 
IMF models seems fully adequate. Moreover, |De Marchi et| 
al. ( 2010 1 conclude that clusters should have born with very 
low values of the MF turnover mass (m c — 0.15 Mq) and it 
was the subsequent dynamical evolution that drove the drift 
of m c to higher values. Since dynamical evolution should 
be especially slow (or virtually absent) in NGC 2419 it is 
likely that a single power-law IMF is a good model for this 
cluster down to very low masses. Indeed, using Eq. 2 by [Del 



Marchi et al. (20101 to obtain rough estimate of m c we get 



m c ~ 0.15 Mq for NGC 2419. 

In the upper panel of Fig. [9] we compare the observed 
LF, in the radial range where it reaches the faintest mag- 
nitudes, with four theoretical LFs from the Dartmouth 
set, with standard values of the physical parameters and 
different exponent of the underlying power-law IMF. For 
Mf814w > 4.5, where the slope of the observed LF appears 
to stabilise, the data points are bracketed between the mod- 
els with MF slope x between -1 and -2. In the following we 
will adopt a x = —1.5 model as the best-fit. 




. Kroupa IMF Padua 
- dH/dM«M- ,J » Basti 
_dN/dM«M" IB Dartmouth 




Figure 9. Observed LF compared with theoretical models, after 
normalisation to the four points in the range 0.5 < Mp%x4,W < 
2.0. Upper panel: comparison with Dartmouth models; different 
line styles are adopted for IMFs with different exponents, as indi- 
cated in the upper left corner of the panel. All the LF models have 
[Fe/H]-2.1, [o/Fe]=+0.4, age=12.0 Gyr. Lower panel: comparison 
with a Basti model with [Fe/H]=-1.8, [ct/Fe]=0.0, age=13.0 Gyr 
and IMF exponent =-1.35 (long dashed line) extrapolated to the 
H-burning limit with a spline (dotted line, see text for details), 
and with a Padua model with [Fe/H]=-1.8, [a/Fe]=0.0, age=13.5 
Gyr and a Kroupa's IMF (dashed-dotted line) . The best fit Dart- 
mouth model is also reported for comparison (continuous line). 



In the lower panel of Fig. [9] we show that a solar-scaled 
composition model with the same global metallicity as the 
standard model, [M/H] = —1.8 (following the prescriptions 
by |Salaris, Chieffi fc Straniero||1993J ), and age=13.0 Gyr 
from the Basti set (dashed line) provides an excellent fit 
to the observed LF over the whole observed range, adopt- 
ing x = —1.35, in good agreement with the best-fit slope 
obtained from the Dartmouth modela 16 ! Basti models are 
limited to m ^ 0.5 Mq, that is fully adequate for the range 
of our LFs but will require some additional assumption to 
extrapolate the total M/L. We took the very simple ap- 
proach to extrapolate the Basti model with a spline (dot- 
ted line) down to the H-burning limit of the best-fit Dart- 
mouth model, with the additional constraint that at that 
point the predicted logA^ must also match. This prescrip- 
tion forces the Basti model to lie below the Dartmouth LF 
for Mfsuw > 7.5, i.e. the regime where stars gives the min- 
imum contribution to the light budget and the maximum 
contribution to the mass budget. For this reason, the M/L 
ratio derived from the integration of the Basti LF should be 
considered as a lower limit. 



13 [M/H] = [Fe/H] + log (0.638 X 1Q["/ Fe l + 0.362), wher e 
[M/H] = \ogZ -logZ G) , 



Salaris, Chieffi & Straniero 



(1993): 



|Ferraro et al.| ( |1999[ l and references therein. With the assump- 
tions listed in Table [I] we obtain [M/H] = -1.8 for NGC 2419, 
corresponding to Z = 0.0003. 

14 http : //193 . 204 . 1 . 62/index . html 

15 http://stev.oapd.inaf . it/cgi-bin/cmd 



16 All the Basti models considered in this paper (either LFs, 
isochrones or integrated properties) have the Reimers' mass loss 
parameter n equal to 0.4, as expected for a cluster with an ex- 
tended blue HB. The choice of this parameter has negligible im- 
pact on our analysis. The Helium abundance of all the considered 
models is Y=0.245, independently of the adopted theoretical set. 
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The overall shape of the observed LF in the range 
2.5 < Mfsuw ^ 4.0 is also acceptably reproduced by 
a Padua model with [M/H] — —1.8, age=13.5 Gyr and 
Kroupa (2001 I IMF. On the other hand it appears to over- 



estimate the counts for Mfsi4W > 4.0, as expected, since 
in this range the Kroupa IMF corresponds to a x = —2.3 
power-law. In this case the choice of an arbitrary power-law 
IMF is not allowed by the web-tool. However we decided 
to also consider this model, as it is expected to provide an 
upper limit to the M/L. The strong peak at Mfsuw — 0.0 
is due to the HB that is included in the Padua LFs and not 
included in those from the other sets of models considered 
here. 

The best-fit Dartmouth model has been plotted also in 
this panel, as a reference. It is interesting to note that all 
the models correctly predict the position of the RGB Bump 



and the star counts in this feature (see Ferraro et al.|1999 
and references therein), at Mfsi4,w — —1.0. It is also worth 
noting that, according to the thorough analysis by Paust et 
al. (20101, a power-law MF with a slope —1.5 > x > —1.0, as 



observed here, is typical of metal-poor GCs like NGC 2419, 
and fully compatible with the range of x covered by clusters 
with similar concentration (see their Fig. 21). 

The three models shown in the lower panel of Fig. [9] 
will be used in the following to estimate the V band 
Mass-to-Light ratio (M/Ly, always expressed in solar units 
[Mq/L vQ ]) of NGC 2419. The Dartmouth LF, that pro- 
vides both a good fit to the observations and a self-consistent 
model down to the H-burning limit will be the reference for 
our best-fit estimate. The other two models presented in the 
lower panel of Fig. [9] will be used for an independent san- 
ity check of our results. They have been explicitly chosen to 
bracket the reference model. Moreover, they differ in details 
of the input physics and theoretical/empirical assumptions. 
Hence they also serve to test the robustness of the derived 
M/Ly to all these factors. 



4.1 The stellar mass-to-light ratio 

The integration of the LF and MF of the best-fit models 
leaves out of the M/L budget two ingredients that may have 
a significant impact on the the final estimate: 

(i) The contribution of the HB and AGB stars, that is 
not included in the Dartmouth and Basti LFs. 

(ii) The contribution of remnants of the stellar evolution, 
i.e. white dwarfs (WD), neutron stars (NS) and black holes 
(BH). 

Concerning the first issue, we have adopted different ap- 
proaches for the different models, again to explore the effects 
of different sources of uncertainty. To include the contribu- 
tion of HB+AGB in the M/L estimate from the Dartmouth 
model we selected all the observed HB and AGB stars using 
the selection box shown in the right panel of Fig. [2] we con- 
verted their magnitudes from F606W to V using the trans- 
formations presented in Fig. IT] and, finally, to solar lumi- 
nosities by assuming M v q =4.83 (seelMaraston & Thomas 
|2000[ ). Then we produced a synthetic HB population with 
the dedicated Dartmouth web toorj roughly reproducing 



the observed morphology, and we derived a mean relation 
to convert V magnitudes into masses. Then summing the 
luminosity and the mass of all the HB+AGB stars we get 
the contributions to be added to the corresponding values 
obtained from the integration of the LF+MF models. It is 
worth to stress here that HB and AGB stars, that may give 
raise to 10 to 20 per cent of a GC V light (depending on age, 
metallicity, IMF and HB morphology of the population), 
provide only a negligible contribution to the mass budget of 
the population. Therefore, the final results are completely 
insensitive to the details of the conversion of their stellar 
magnitudes into masses. To make the HB+AGB contribu- 
tion consistent with the integration of the MS+RGB LF 
models we (a) normalise the LF models to the observed LF 
within the half-light radius before performing the LF/MF 
numerical integration, and (b) we only consider AGB and 
HB stars lying in the same region. This gives the additional 
advantage that multiplying the derived luminosity by two 
we can obtain an estimate of the total cluster luminosity, 
and consequently My, that can be compared with the fully 
independent estimates available in the literature (see Tablefl] 
and |Belrazzim|2007| ). 

For the Basti model we proceeded as follows. We pro- 
duced a synthetic CMD with the dedicated web tool, for 
a simple stellar populations (SSP^rj with age=13 Gyr, 
[M/H]=-1.8, solar scaled composition and power-law IMF 
with x = —1.35. We use the synthetic population to com- 
pute the fraction of the total V light and mass contributed 
by stars in the HB and AGB phases. It turns out that these 
stars provide 18.2 per cent of the total V light and 0.1 per 
cent of the total stellar mass. The synthetic population has 
indeed a Blue HB morphology but not so extreme as the 
observed one. For instance, the faintest (bluest) HB star of 
the synthetic population has Mi ~ 1.5, while the observed 
HB distribution of NGC 2419 extends down to Mi ~ 5.5. 
Hence the derived fraction should be considered as an upper 
limit to the light contribution of these stars, that is appro- 
priate since, in the present context, the Basti model has been 
chosen to provide a lower limit to M/Lv- 

Finally, it should be noted that for the Padua models 
the implicit assumption is that all the HB stars lie in a red 
clump, very different from the actual HB morphology of the 
cluster. However also this approximation should have only 
a minor impact on the final M/Ly ratio since, as said, the 
overall contribution of HB stars to the light budget is ^ 20 
per cent, hence, once the number of HB stars is fixed, vari- 
ations in their luminosity distribution should induce varia- 
tions in the overall contribution significantly smaller than 
this, i.e. of order of a few per cent. 

Turning to point (ii), the remnants contribute only to 
the mass budget, since even WDs (the only class of non-dark 
remnants) are so faint that their contribution to the total 
light is negligible. Basti models of the integrated proper- 
ties of SSPs ( Percival et al.|2009 l provide the mass fraction 
in remnants obtained by integrating a Kroupa ( 2001 1 IMF 



http://stellar.dartmouth.edu/ models/shb.html 



with appropriate mass limits and adopting simple initial- 
final mass relations. For a SSP with age= 13 Gyr and 



18 A simple stellar population is an ensemble of stars sharing the 
same age, IMF and chemical composition |[Renzini fc Fusi Pecci| 



1988) 
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[M/H] = —1.8, the fraction of the total cluster stellar mass 
in BH, NS, and WD is 0.03, 0.008, and 0.139, respectively. 
All the remnants are produced by processes involving explo- 
sions that are generally believed to impart a velocity kick 
that may lead to the ejection of the object from the cluster 
(see, e.g. Pfahl et al. 2002). In a recent analysis Kruijssen 



M/Lv = 1.5, with a typical 1-er error ~ 0.1, accounting for 
observational and theoretical uncertainties. 

The M/Lv values reported above are lower than the 



(20091 estimates that globular clusters retain 0.3-7 per cent 



of the BH and 0.1-4 per cent of the NS. Combining the mass 
fractions with the retention factors it should be concluded 
that the actual contribution of BH and NS to the present- 
day mass budget of the cluster is ^0.2 per cent and ^0.03 
per cent, respectively. This is clearly negligible and we do 
not consider it in our computation of M/Lv- For WD we 



adopt (conservatively) a retention factor of 0.75 (from Krui- 
jssen| [2009), so obtaining a total contribution of WD to the 
present-day mass budget of the cluster of 10 per cent. 

Given all the above assumptions, integrating the theo- 
retical V band LF/MF of the Dartmouth model that best 
fits the observed LF (Fig. [9} from the RGB Tip down to the 
Hydrogen-burning limit, and including the contribution of 
HB+AGB stars and dark remnants as described above, we 
obtain M/Lv =1.49. The associated value of the absolute 
integrated V magnitude is My = —9.6, in good agreement 
with the independent estimates based of surface photom- 
etry (M v = -9.5, Table [l). Adopting models with IMF 
exponent x = —1.0 or x = —2.0, instead of the best-fit 
value x = -1.5, we obtain M/L v =1.11, and M/L v =2.09, 
respectively. This gives an idea of the maximum effect on 
M/Lv that can be obtained by changing the IMF slope, 
since x values outside the —2.0 < x < —1.0 range are 
clearly excluded by the observed LF (at least in the acces- 
sible range of magnitudes) . It is interesting to note that the 
assumption of x = —1.0 or x = —2.0 imply total luminosities 
that are hardly compatible with integrated photometry, i.e. 
Mv = —10.1, and Mv = —9.0, respectively. Adopting a 100 
per cent retention fraction of all the dark remnants, the final 
ratio moves only from M/Lv =1.49 to M/Lv =1.60. As a 
further consistency check we integrated the I band LF/MF 
of the same model, then converting the total I luminosity 
into V luminosity by adopting (V — I) q = 0.688, from 
Holmberg, Flynn fc Portinari| ( |2006| , and' (V - I) int = 1.05 
as the observed integrated color of NGC 2419, from |Peter-| 

( 1993 1. With this approach we find M/Lv =1.55 for the 



best-fit model, fully consistent with the result obtained from 
the V integration. Since the normalization to the observed 
LF within tk has been performed in the I (or equivalently 
F814W) band, this path of integration is expected to pro- 
vide the most reliable estimate of the total luminosity. We 
find Mv = —9.5, in excellent agreement with the estimate 
from integrated photometry. 

The direct V(I) integration of the Basti and Padua 
models, with the prescriptions described above for remnants 
and HB+AGB stars, gives V-band mass-to-light ratios of 
M/Lv=l. 18(1.24) and M/Z,y=1.59(1.72), respectively. The 
derived values bracket the best-fit estimate from the ex- 
pected sides (i.e. Basti= lower limit, Padua= upper limit). 
The absolute integrated magnitude from the integration in 
I gives M v = -9.6 (Basti) and M v = -9.8 (Padua), both 
compatible with integrated photometry. 

We conclude that the total stellar M/Lv ratio should 
lie in the range 1.2 < M/Lv < 1-7, with a best-fit value 



independent estimates obtained by Baumgardt et al. ( 2009 
M/Lv = 2.05) and |Ibata et al.] ( |2011a|b[ M /L v = 1.90) 
from the kinematics of cluster stars (see also McLaughlin 
|fc van der Marel"||2005| and references therein), but clearly 
not incompatible with therrpj We defer the detailed dis- 
cussion of a possible tension between the stellar and the 
dynamical M/L to a future contribution, dedicated to the 
analysis of dynamical models of the cluster including dark 
matter (Nipoti et al., in preparation). On the other hand, 
the agreement with predictions from SSP models is satisfac- 
tory ( |Ibata et aL]|2011a l. 

The above results are strongly dependent on the as- 
sumption that the slope of the MF from our last measured 
point at m ~ 0.5 Mq (Mfsi4w = 6.4) to the H-burning 
limit, at m ~ 0.1 Mq (Mfsi4w — 11.0) does not vary. 
While it is quite likely that a single power-law is a good rep- 
resentation of the cluster MF down to very low masses, it is 
worth exploring the effects of a turnover at low masses ( |De| 
Marchi et aT|2010[ [Bastian, Covey fc Meyer||2010[ ). With a 
very conservative approach, given the discussion of Sect. HI 
we consider the (unrealistically extreme) case that the clus- 
ter does not contain stars with mass lower than 0.3 Mq. Un- 
der this condition, we obtain M/Lv = 0.89, 0.85, and 1.14 
from the Dartmouth, Basti and Padua models, respectively. 
Thus, even adopting the most conservative assumptions, a 
robust lower limit to the stellar mass to light M/Lv > 0.8 
is obtained. 

As a final consideration it may be worth noting that the 
above estimate neglects the presence of binaries, while it is 
likely that a significant fraction of these systems should be 
there (likely ^ 20 per cent, by analogy with other GCs with 
similar density, see |Sollima et"aT1|2007| |Milone et aLl|2012| 
Ibata et al.||2011a and discussion therein). At the distance 
of Galactic globulars, all binaries are unresolved and are 
necessarily included in LFs as single stars slightly brighter 
than the primary of the system (up to — 2.51og2 ~ —0.75 
mag, in case of equal mass components). This leads to two 
effects, both concurring to obtain M/L estimates that are 
lower than the true values: (a) low-mass stars move from 
their natural magnitude bin to the brighter bin including 
the brightened primary, thus leading to an artificial flatten- 
ing of the M F ( |Malkov fc Zinnecker|2000| |Bastian, Covey fc 
Meyer|20~10 1 , and (b) the unresolved binary appears to have 
a M/L that is lower than the M/L of each of the individual 
components. Hence, neglecting binaries in our analysis does 
not have any effect on the derived lower limits on M/Lv- 
Some simple simulation^] indicate that the true M/Lv can 



19 It should also be considered that the uncertainties in the dis- 
tance and reddening estimate are not included in our error budget. 
However these are not expected to have a major impact on our es- 
timates, since the proper match between observed and predicted 
LFs is also constrained by distance-independent and reddening- 
independent features, like the difference in magnitude between 
the RGB bump and the rapid rise of star counts before the TO 
point. 

20 Performed with synthetic populations including binaries ex- 
tracted from a flat distribution of mass ratios ||Fisher et al~ 2005), 
following the same approach as in Ibata et al. 1(201 la I. 
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be larger than our best estimate by a factor from ~ 1.05 to 
~ 1.25 for binary fractions going from 10 per cent to 50 per 
cent, i.e. an amount comparable with the overall uncertainty 
in the M/Ly estimate. This effect should help to reconcile 
the derived stellar M/Ly values with the dynamical esti- 
mates. For example, assuming a binary fraction of 20 per 
cent, M/Ly is expected to grow from ~ 1.5 to ~ 1.7. 



5 SUMMARY AND CONCLUSIONS 

We have used archival HST-WFC3 data to obtain the deep- 
est Luminosity Function so far for the massive and remote 
globular cluster NGC 2419. The derived LF was corrected 
for the effects of incompleteness (and its strong radial vari- 
ation) and bin migration. Comparing the LFs obtained in 
different radial ranges, also including the wide-field ground- 
based dataset by D08, we find that the shape of the LF is 
the same all over the cluster, at least in the range of lumi- 
nosity covered by our data. This result strongly supports the 
conclusion by D08 that there is no mass segregation within 
the cluster, because of the inefficiency of two-body relax- 
ation in establishing the condition of energy equipartition 
(see Baumgardt et al.|2009 for a detailed discussion). As a 
consequence, this supports the validity of the hypothesis of 
constant stellar mass-to-light ratio adopted by [Baumgardt 
et al. (20091 and Ibata et al. (2011a b) in their analysis of 
the cluster kinematics. 

We have fitted our observed LF with theoretical models 
from three different sets. In the range covered by our data, a 
power-law IMF with exponent x ~ —1.5 provides a reason- 
able fit to observations, in agreement with the trends found 



not seem to support this scenario ( Bastian, Covey & Meyer 



20101 



by |Paust et aL] ( |2010[ ). Integrating the LF/MF of the best fit 
models and extrapolating the results down to the H-burning 
limit, with suitable assumptions to account for the contri- 
bution of HB+AGB stars and stellar remnants, we find a V 
band mass-to-light ratio in the range 1.2 < M/Ly < 1.7, 
with a best-fit estimate M/Ly = 1.5 ± 0.1, slightly lower 



than the most recent dynamical estimates (Baumgardt et al. 
Ibata et al.|2011a|bl, but still compatible with them, 



2009 



when all the sources of uncertainty are taken into account 
On the other hand, assuming that the cluster has no star 
with m ^ 0.3 Mq, we obtain a robust and conservative 
lower limit M/Ly > 0.8. 

In this context it is worth mentioning two astrophysi- 
cal processes that may provide a way to increase the mass 
fraction of dark remnants, thus helping to reconcile the stel- 
lar and the dynamical M/L estimates. First, all the avail- 
able models for the formation of multiple stellar populations 
in GCs invoke progenitors with total masses much larger 
than present-day clusters ( D'Ercole et al.|2008 Decressin et 



[aT||2010| |Bekki|[20TT| ). In particular, within the theoretical 
framework they adopt to explain the claimed large spread 
of He abundance in NGC 2419, Idi Criscienzo et al.l (|2011b I 



conclude that the mass of the cluster at its birth should 
have been from 15 to 30 times the present-day mass. The 
resulting deep potential well in the first 50-100 Myr of the 
cluster life may have helped in retaining a larger fraction of 
dark remnants within the inner regions of the proto-cluster. 
Second, a top-heavy IMF ( Marks et al.|2012 | may enhance 
the contribution of dark remnants to the mass budget of the 
cluster, but currently available observational evidence does 



The fact that mass-to-light ratio derived from SSP mod- 
els are typically larger than those found from dynamical 
analysis is likely accounted for by the selective loss of low- 



mass stars due to the dynamical evolution of clusters ( Krui- 
jssen & Mieskc 20091. The present study suggests that the 
comparison can be pushed to an higher level, since, on one 
hand dynamical estimates not simply based on the central 
velocity dispersion are becoming available (see, e.g. Lane 



et al. 20111, and on the other hand there is the observa- 
tional material ( |Paust et al.|2010| |De Marchi et al.|20l6"| ) to 
derive very robust stellar mass-to-light ratios by integrat- 
ing observed LF/MF nearly reaching the H-burning limit. 
These should reflect the present-day stellar content of clus- 
ters, not just an educated guess based on average properties 
of idealized SSP models. These measures based on actual 
star counts would provide much more stringent (and less 
model-dependent) constraints on the stellar M/L ratio, to 
be compared with dynamical estimates. 
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